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Abstract. We develop a new variational formulation of the inverse Stefan problem, where information on the 
heat flux on the fixed boundary is missing and must be found along with the temperature and free boundary. 
We employ optimal control framework, where boundary heat flux and free boundary are components of the 
control vector, and optimality criteria consist of the minimization of the sum of L 2 -norm declinations from 
the available measurement of the temperature flux on the fixed boundary and available information on the 
phase transition temperature on the free boundary. This approach allows one to tackle situations when the 
phase transition temperature is not known explicitly, and is available through measurement with possible 
error. It also allows for the development of iterative numerical methods of least computational cost due to 
the fact that for every given control vector, the parabolic PDE is solved in a fixed region instead of full 
free boundary problem. In Inverse Problems and Imaging, 7, 2(2013), 307-340 we proved well-posedness in 
Sobolev spaces framework and convergence of time-discretized optimal control problems. In this paper we 
perform full discretization and prove convergence of the discrete optimal control problems to the original 
problem both with respect to cost functional and control. 

Key words: Inverse Stefan problem, optimal control, second order parabolic PDE, Sobolev spaces, 
energy estimate, embedding theorems, traces of Sobolev functions, method of finite differences, discrete 
optimal control problem, convergence in functional, convergence in control. 

AMS subject classifications: 35R30, 35R35, 35K20, 35Q93, 65M06, 65M12, 65M32, 65N21. 


1 


1 Description of Main Results 
1.1 Introduction and Motivation 

Consider the general one-phase Stefan problem ([151 126]): hnd the temperatnre fnnction 


u{x,t) and the free bonndary x = s{t) from the following conditions 

{a{x,t)ux)x + b{x,t)ux + c{x,t)u — Ut = f{x,t), for {x,t)Efl (1.1) 

u{x, 0) = (j){x), 0 < X < s(0) = So (1.2) 

a(0, t)ux{0, t) = g{t), 0 <t <T (1.3) 

a{s{t),t)ux{s{t),t) + 7 (s(t), t)s'(t) = x{s{t),t), 0<t<T (1.4) 

u{s{t),t) = 0<t<T (1.5) 

where a, b, c, f, 0 , g, 7 , X; h known fnnctions and 

a{x,t) > oo > 0, So > 0 (1.6) 

12 = {{x,t) : 0 < X < s{t), 0 < t < T} 


In the physical context, / characterizes the density of the sonrces, 0 is the initial temperatnre, 
g is the heat flux on the hxed boundary and /i is the phase transition temperature. 

Assume now that some of the data is not available, or involves some measurement error. 
For example, assume that the heat flux g{t) on the hxed boundary a: = 0 is not known and 
must be found along with the temperature u{x,t) and the free boundary s{t). In order to 
do that, some additional information is needed. Assume that this additional information is 
given in the form of the temperature measurement along the boundary a: = 0: 

u{0,t) = for 0 < t < T (1.7) 

Inverse Stefan Problem (ISP): Find the functions u{x,t) and s{t) and the boundary heat 
flux git) satisfying conditions / li.il) - ((i. ?p . 

Motivation for this type of inverse problem arose, in particular, in the modeling of bio¬ 
engineering problems on the laser ablation of biological tissues through Stefan problem fll.lD - 
fll.6p . where s{t) is the ablation depth at the moment t. The boundary temperature mea¬ 
surement u{0,t) contains an error, which makes it impossible to get reliable measurement of 
the boundary heat hux g{t), and the ISP must be solved for its identihcation. This approach 
allows us to regularize an error contained in a measurement i^(t). Another advantage of this 
approach is that, in fact, condition fll.5l) can be treated as a measurement of the temperature 
on the ablation front, and our approach allows us to regularize an error contained in tem¬ 
perature measurement fi{t) on the ablation front. Still another important motivation arises 
in optimal control of the Stefan problem, where controlling g{t) is equivalent of controlling 
external temperature along the hxed boundary. It should be pointed out that the method 
of this paper can be applied to different type of inverse problems. For example, fll.7p can be 
replaced with 

u{x,T) = w{x), for 0 < a; < s(T), 

meaning that measurements are taken for the hnal temperature distribution w{x) and dual 
ablation depth s(T). Instead of identihcation of the boundary hux g, one can consider the 
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inverse free boundary problem with any of the unknown coefficients a, b, c or right hand side 

/• 

The ISP is not well posed in the sense of Hadamard. If there is no coordination between 
the input data, the exact solution may not exist. Even if it exists, it might be not unique, 
and most importantly, there is no continuous dependence of the solution on the data. The 
ISP was first mentioned in ra. in the form of Ending a heat flux on the fixed boundary 
which provides a desired free boundary. This problem is similar to a non-characteristic 
Cauchy problem for the heat equation. The variational approach for solving this ill-posed 
inverse Stefan problem was performed in |7118]. The first result on the optimal control of 
the Stefan problem appeared in [36] . It consists of Ending the optimal value of the external 
temperature along the fixed boundary, in order to ensure that the solutions of the Stefan 
problem are close to the measurements taken at the final moment. In [36], the existence 
result was proved. In [5B] the Frechet differentiability and the convergence of the difference 
schemes was proved for the same problem and Tikhonov regularization was suggested. Later 
development of the inverse Stefan problem was along these two lines: Inverse Stefan problems 
with given phase boundaries were considered in [21IH El El IHl 1121 CSl ESI E2l ESj; optimal 
control of Stefan problems, or equivalently inverse problems with unknown phase boundaries 
were investigated in [H [191 1201 1211 1221 1231 ESI 1291 (23 EQI [SH [35l [16]. We refer to 
monography [16] for a complete list of references of both types of inverse Stefan problems, 
both for linear and quasilinear parabolic equations. The main methods used to solve the 
inverse Stefan problem are based on variational formulation, method of quasi-solutions or 
Tikhonov regularization which takes into account ill-posedness in terms of the dependence 
of the solution on the inaccuracy involved in the measurement fll.7p , Frechet differentiability 
and iterative conjugate gradient methods for numerical solution. Despite its effectiveness, 
this approach has some deficiencies in many practical applications: 

• Solution of the inverse Stefan problem is not continuously dependent on the phase 

transition temperature small perturbation of the phase transition temperature 

may imply signiflcant change of the solution to the inverse Stefan problem. Accordingly, 
any regularization which equally takes into account instability with respect to both z/(t) 
from measurement fll.7p . and the phase transition temperature /i(t) from fll.51) will be 
preferred. It should be also mentioned that in many applications the phase transition 
temperature is not known explicitly. In many processes the melting temperature of 
pure material at a given external action depends on the process evolution. For example, 
gallium (Ga, atomic number 31) may remain in the liquid phase at temperatures well 
below its mean melting temperature ([26]). 

• Numerical implementation of iterative gradient type methods within the existing ap¬ 
proach requires solving the full free boundary problem at every step of the iteration, 
and accordingly has quite a high computational cost. An iterative gradient method 
which requires solution of the boundary value problem in a fixed region at every step 
would definitely be much more effective in terms of the computational cost. 

The main goal of this project is to develop a new variational approach based on the 
optimal control theory which is capable of addressing both of the mentioned issues and 
allows the inverse Stefan problem to be solved numerically with least computational cost 
by using conjugate gradient methods in Hilbert spaces. In [T] we proved the existence 
of the optimal control and convergence of the family of time-discretized optimal control 
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problems to the continuous problem. In this paper we perform full discretization through 
hnite differences and prove the convergence of the discrete optimal control problems to the 
continuous problem both with respect to cost functional and control. We employ Sobolev 
spaces framework which allows us to reduce the regularity and structural requirements on 
the data. We address the problems of Frechet differentiability and application of iterative 
conjugate gradient methods in Hilbert spaces in an upcoming paper. 

Throughout the paper we use the usual notation for Sobolev spaces according to references 
[2n ISl [28l |33l 131]. Notation is described below in Section 11.21 


1.2 Notation of Sobolev Spaces 

L2[0,T] - Hilbert space with scalar product 


(m,u) 



uvdt 


W2[0, T],k = 1, 2,... - Hilbert space of all elements of T2[0, T] whose weak derivatives up to 
order k belongs to L 2 [ 0 ,T] and scalar product is dehned as 


(u,u) = 

Jo s=0 


d^u d'^v 
dd dt^ 


dt 


W2‘[0,T] - Banach space of all elements of L2[0,T] with hnite norm 


'“Vyo,Ti^ + 


/ dt / 

'0 Jo 


|M(t) -u(r)p 


\t-T 


■dr 


L 2 (H) - Hilbert space with scalar product 


(u, v)= uvdxdt 

Jq 

W 2 ^’°(H) - Hilbert space of all elements of L 2 {kl) whose weak derivative belongs to ^2(11), 
and scalar product is dehned as 

- Hilbert space of all elements of L 2 (H) whose weak derivatives ^ belong to 
L 2 (H), and scalar product is dehned as 

du dv du dv' 


V 2 (f^) - Banach space of all elements of hF 2 ^’°(H) with hnite norm 


Q \ OX ox ot ot / 


u\\v 2{n) = (esssupo<t<T\\u{x,t)\\L^[o,sit)] + 


du 


dx 


L2{n) 
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V"2^’°(f2) - Banach space which is the completion of in the norm of V 2 (ff). It consists 

of all elements of V 2 (fl), continnous with respect to t in norm of L 2 [ 0 ,s(f)] and with hnite 
norm 



( max 
V 0 <t<T 


'“(^)^)llL[0,s(i)] + 


du 

dx 



1 

2 


W 2 '^{yL) - Hilbert space of all elements of -^ 2 (^ 1 ) whose weak derivatives |^, belong 

to 1/2 (H), and scalar product is dehned as 


(n,n) 



du dv du dv 
^ dxdx^ dt dt 


di^u di^v 

dx‘i Qx‘^ 


dxdt 


1.3 Optimal Control Problem 

Consider a minimization of the cost functional 


J{v) =/3o||M(0,f) - z/(t)||i 2 [ 0 ,T] + -hWllLp.r] (1-8) 

on the control set 


Vr = {v = {s,g) e W^[0,T] X W^[0,T] : 5 < s{t) < l,s{0) = So,s'{0) = 0, 

max( ||s||^ 2 ; ||^||^i < R} 

where 6,1, R, (3o, (3i are given positive numbers, and u = u{x,t;v) be a solution of the Neu¬ 
mann problem fll.ip - fll.4l) . 

Definition 1.1 The function u G Wl’^iVL) is called a weak solution of the problem U.l\} - 
^1-4^ ifu{x,0) = 0(x) G iy 2 ^[ 0 ,So] and 

j>T i>s{t) 

0 = / / [aUx^x — hux^ — CM<h -|- Ut^ -|- /<h] dx dt 

Jo Jo 

+ f b{s{t),t)s'{t)-x{s{t),t)]^{s{t),t)dt+ [ g{t)^{0,t)dt (1.9) 

Jo Jo 

for arbitrary <h G 


We also need a notion of weak solution from V 2 (f^) of the Neumann problem: 
Definition 1.2 The function u G V 2 (f^) is called a weak solution of U.l\)-^1.4\) if 


rT psit) 


'■so 


0 = 
rT 


[aux^x — bux^ — cn<h — u^t + f^] dxdt — / (j){x)^{x, 0) dx+ 


'0 JO 


g{t)^{0,t) dt+ / [y{s{t),t)s'{t)-u{s{t),t)s'{t)-x{s{t),t)]^{s{t),t)dt (1.10) 


>0 


for arbitrary <h G such that = 0. 

If M is a weak solution either from V 2 (f^) (or kh2^’^(f2)), then traces u\^^q and are 

elements of L 2 [ 0 ,T], when s G kh 2 ^[ 0 ,T] ([2Hll2l|) and cost functional J7’(n) is well dehned. 
Furthermore, formulated optimal control problem will be called Problem I. 
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1.4 Discrete Optimal Control Problem 


Let 

= {tj =j-T, j = 0,1,..., n} 

be a grid on [0,T] and t = ^. Consider a discretized control set 

K = bin) e : 0 < i < s, < /, mfK(||[s]„|£,; |||!,],.|£.) < R^] 

where, 


L,2 — 


[-S]n — (-So, -Si, ..., Sn) G [g]n — {go, Ql, •••, Qn) G 

n—1 n n—1 n—1 n 


k- 


k=0 k=l k=0 k=0 k=l 

where we assign s_i = Sq and use the standard notation for the hnite differences: 

_ -Sfc -Sfc—1 ^ _ -Sfc+i -Sfc ^ _ -Sfc+1 2sfc -|- -Sfc—i 

-Slfc ~ ) -St,fc = , -Sit fc = 5 • 

r T 

Introduce two mappings Qn and Vn between continuous and discrete control sets: 

Qn{v) = [v]n = ([-s]n, [fi']n), for V E Vr 
where Sk = s(tk),gk = g(tk), k = 0,1,..., n. 

Pni[v]n) =V^= (s^ g^) E W^[0, T] X T] for [v]n E V^, 

where 

^ I So + 0 < t < r, _ 

\ Sk-i + (t - 4-1 - ^)sik-i + |(^ - ik-ifs^^k-i 4-1 <t<tk,k = 2,n. 

g'^it) = fl'fc-i + ——- 4-i), 4-1 <t<tk,k = 1^. 


( 1 . 11 ) 


Let us now introduce spatial grid. Given [t;]„ E V^, let {po,Pi, ■ ■ ■ ,Pn) be a permutation of 
(0,1, • • • ,n) according to order 


Spo — Spj ^ ^ -Sp„ 

In particular, according to this permutation for arbitrary k there exists a unique jk such that 

Sk = Sp.^ (1.12) 

Furthermore, unless it is necessary in the context, we are going to write simply j instead of 
subscript jk. Let 

(.^po = {xi : Xi = i ■ h, i = 0,1,..., 

be a grid on [0, Sp^] and h = Furthermore we always assume that 

niQ 

h = 0(\/r), as r ^ 0. (1.13) 

We continue construction of the spatial grid by induction. Having constructed cUp^. -^ on 
[0, -Spj,_ J we construct 

Wpfc = {xi-. i = 0,1,-- - ,mf^} 
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on [0, Spf,], where and ineqnality is strict if and only if Sp^ > Sp^_^] for i < 

points Xi are the same as in grid 0 Jpi^_^. Finally, if Sp^ < I, then we introduce a grid on [sp„, /] 

U = {xi : Xi = Sp^ + {i - m^n^)h, i = ■■ ■ ,N] 


of stepsize order h, i.e. h = 0{h) as h ^ 0. Furthermore we simplify the notation and write 
= rrifc. Let 

hi = Xi+i - Xi, i = 0,1, - ■ ■ ,N - 1] 

and assume that 


ruk —)■ +C) 0 , as n ^ cxD. 


Introduce Steklov averages 


dk{x) = 




T 


d{x, t) dt, hk = 


f*^k 


h{t) dt, dik = 




tk-l 


'ifc-1 


hiT 


d{x, t) dt dx. 


^k — 1 


where z = 0,1, • • • , iV — 1; k = 1, ■■■ ,n-, d stands for any of the functions a, b, c, f, and 
h stands for any of the functions u, /z, g or g^. Given v = {s,g) G Vr we dehne Steklov 
averages of traces 

Xs = -[ X{s{t),t)dt, {jss')’'= - [ j{s{t),t)s'{t)dt. (1.14) 


Given [n]„ = ([s]n, [g]n) £ Vr we define Steklov averages and (7s"(s”)')^ through fll.ldp 
with s replaced by s” from fll.lip . 

Let 0” be a piecewise constant approximation of (j): 


0"'(x) = (pi := 0(xj), for Xi < X < Xj+i, z = 0,iV — 1 
Next we define a discrete state vector through discretization of the integral identity fll.Qp 


Definition 1.3 Given discrete control vector [v]n, the vector function 

[u{[v]n)]n = (zz(0), zz(l),zz(n)), zz(fc) G A; = 0, • • • , n 

is called a discrete state vector if 

(a) First tuq + 1 components of the vector zz(0) G satisfy 

Ui{0) = pi := p{xi), z = 0,1, • • • , mo; 


(b) Recalling fll.l2p . for arbitrary k = !,■■■ ,n first mj + 1 components of the vector 
u{k) G solve the following system of mj + 1 linear algebraic eguations: 


hh 


flofc + hbok — h Cofc H- 

r 


Uo{k) 


o-ok + hb{ 


Ok 


h^ 


ui{k) = —uo{k - 1) - hVofc - hg^, 

T 

r h^h- i' 

^i—\^khiUi—\{fp “ 1 “ ^i—l^khi aikhi_-i bikhihi_-i Cikhihi_Y -\- 


T 


Ui{k)- 




hph- 

Ui+i{k) = -h^ihi_ifik + ' - 1), z = 1, 


r 


,mj — 1 


~0‘mj — l,kUmj — l{ki) 4“ Um, —(A^) hmj — l (7s” ) ) X. 


(1.15) 
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(c) For arbitrary k = 0,1,...,n, the remaining components of u{k) G are calculated 


as 


Ui{k) = u{xi] k), mj < i < N 

where u{x; k) G 11 ^2^[0, /] is a piecewise linear interpolation of {ui{k) : i = 0, - , mj}, 

that is to say 

7\ n \ \ ~ \ ^ ^ • n 1 

u[x\ k) = Ui{k) H--- [X — Xij, Xi < X < Xi+i,t = 0, • • • ,mj — 1, 

fli 


iteratively continued to [ 0 , /] as 


r/l 


u{x; k) = uifl'^Sk — x;k), 2 ” '^Sk < x < 2 "'Sfc, n = 1, n^, < n* = 1 + log 2 7 (1-16) 


L5. 


where [r] means integer part of the real number r. 


It should be mentioned that for any k = 1,2, •• • ,n system (I3.6ip is equivalent to the 
following summation identity 


ra-i — l 


^ ^ fik^i 


i=0 


+ 




Vmj + QkVo = 0, 


(1.17) 


for arbitrary numbers pi,i = 0,1, ■ ■ ■ , mj. 

Consider a discrete optimal control problem of minimization of the cost functional 

n ^ n 2 

'In{[v]n) = {uo{k) - {k) - fik) (1-18) 


k=l 


k=l 


on a set Vf^ subject to the state vector dehned in Dehnition 1.3. Furthermore, formulated 
discrete optimal control problem will be called Problem 

Throughout, we use piecewise constant and piecewise linear interpolations of the discrete 
state vector: given discrete state vector [u([u]„)]„ = (m( 0 ),-u( 1 ), ...,u{n)), let 


u'^{x, t) = u{x; k), if tk-i < t < tk, 0 < x < I, k = 0,n, 


iF{x, t) = ufx] k — 1) + uj{x-, k){t — tk-i), if 4 -i < t <tk, 0 < x < I, k = l,n, 
iF{x, t) = u{x-, n), if f > T, 0 < x < /. 


u'^{x, t) = Ui{k), if tk-i < t < tk, Xi < X < Xj+i, k = l,n, i = 0, N — 1. 
Obviously, we have 

u" G 14(T>), h" G WI’\D), G L2{D). 

As before, we employ standard notations for difference quotients of the discrete state vector: 

Ui+i{k) - Ui{k) Ui{k) - Ui{k - 1) 

Uix[k) = ---, Uii = -, etc. 


hi 


T 


















1.5 Formulation of the Main Result 

Let 

D = {{x,t) : 0 < X < I, 0 < t < T} 

Throughout the whole paper we assume the following conditions are satisfied by the data: 

a,b,c ^ Loo(D), f G L 2 {D), 

0 e hL2^[0, So], 7) X ^ /i, z/G L2[0, T], 

the coefficient a satisfies (ESI) almost everywhere on D, the weak derivatives || exists 
and 

dt < +CX0. (1-19) 

Our main theorems read: 

Theorem 1.1 The Problem I has a solution, i.e. 

V, = {veVR- J{v) = X= inf J{v)} ^ 0 

vGVr 

Note that Theorem 11.11 was already proved in jl] by using method of lines. 

Theorem 1.2 Sequence of discrete optimal control problems approximates the optimal 
control problem I with respect to functional, i.e. 

lim In, = J*, (1.20) 

n^+oo 

where 

In, ^ 1)2, ... 

If G Vfl is chosen such that 

h ^n([^]ne) ^ '^rit T 0, 

then the sequence Vn = {sn,gn) = 'Pn{[v]nJ converges to some element n* = (s*, 5 f*) G 14 
weakly in hT2^[0, T] x hT2^[0, T], and strongly in 1T2^[0, T] x L2[0, Tj. In particular Sn converges 
to s* uniformly on [0,T]. For any 5 > 0, define 

= 12* n {x < s* (t) — 5, 0 < t < T} 

Then the piecewise linear interpolation vP of the discrete state vector [^[njnjn converges to 
the solution u{x,t-,v^) G of the Neumann problem U.lp - Ihl.Rl ) weakly in 

2 Preliminary Results 

In Lemma ITT] below we prove existence and uniqueness of the discrete state vector [n([n]n)]n 
(see Definition 11.31) for arbitrary discrete control vector [n]„ G Vff. In Lemma [2.21 we remind 
a general approximation criteria for the optimal control problems from ([37]). In Lemma lT3] 
we recall some properties of the mappings Qn and Vn between continuous and discrete control 
sets. 


da r ^ r 


CSSSUP()-^x<l 
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Lemma 2.1 For sufficiently small time step t, there exists a unique discrete state vector 
[M([w]n)]n for arbitrary discrete control vector [v]n G Vf. 

Proof. As it is mentioned above, for any k = 1,2, •• • ,n system fl3.6ip is equivalent to 
the summation identity fll.lZp for arbitrary numbers r]i,i = 0,1, • • • ,mj. Let {ui{k)} be a 
solution of the homogeneous system related to 03.611) . i.e. 

3t = (7.»(s")')‘ = = U = Mk - 1) = 0. 

By choosing in 01.171) pi = Ui{k) we have 


m-i — l 


rrin—l 


hiaikul^{k) + - (^) = hikUi^{k)ui{k) + Cikul{k) ( 2 . 1 ) 


i=0 


i=0 


j=0 


Using 01.6p and Cauchy inequality with e > 0 we derive that 


m-,' — 1 


m-,' — 1 


m-i — l 


ao ^ hiuf{k) + ^ ^ hiul{k) < ^ ^ ML(^) + (^ + ^) (2-2) 


r — 2 

2=0 2 = 0 2=0 


2=0 


where 


M = max (||a||L^(ij); ||6 ||l^(d); ||c||l^(d)). 
By choosing e = oq/M in 02.21) we have 


rrij —1 

y ^ +(““—) hiU^ik) < 0, 

i=o ^ 


i=0 


(2.3) 


where 

) ' 

From 02.3p it follows that Ui{k) = 0, i = 0,1, • • • ,mj, and hence the homogeneous system 
only has a trivial solution for r < tq. Accordingly, system is uniquely solvable and therefore, 
for any given discrete control vector [v]n there exists a unique discrete state vector defined 
by Dehnition 1.3. Lemma is proved. 

The following known criteria will be used in the proof of Theorem 11.21 


Lemma 2.2 Sequence of discrete optimal control problems In approximates the contin¬ 
uous optimal control problem I if and only if the following conditions are satisfied: 

(1) for arbitrary sufficiently small e > 0 there exists number Ni = Nfie) such that Qn{v) G 
Vfi for all V G Vr-^ and N > Ni; and for any fixed e > 0 and for all v G VR-e the 
following inequality is satisfied: 


limsup (Zn{Qn{v)) - J{v)) < 0. (2.4) 

N^oo y ' 


(2) for arbitrary sufficiently small e > 0 there exists number N 2 = N 2 {e) such thatVN{[v]N) G 
Vr+e for all [u]Ar G and N > N 2 ; and for all [uJat G , N > 1 the following in¬ 
equality is satisfied: 

\imsup(j{VNi[v]N))-ZN{[v]N)) <0. (2.5) 

N^oo k / 
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( 2 . 6 ) 


(3) the following inequalities are satisfied: 

limsupj7*(e) > liminf J7*(—e) < 

where jA^e) = inf Jiu). 

Vr±, 

Next lemma demonstrates that the mappings Qn and Vn introduced in Section 11.41 satisfy 
the conditions of Lemma 12.21 

Lemma 2.3 (f^) For arbitrary sufficiently small e > 0 there exists such that 

Qn{v) G Vfi, for all v G Vr-e and n > n^. (2.7) 

Vni[v]n) G Vr+e, for all [v]n G Vfi and n > n^. (2.8) 

Proof. Let 0<e<<i?, nG Vr-^ and Q{v) = [v]n = ([s]n, [9\n)- By applying Cauchy- 

Bunyakovski-Schwarz (CBS) inequality and Fubini’s theorem we have 


n—1 


fc=i 


n-l ^ *'=+1 


E 

k=l 


t3 l 


S'{t) 




s'{t 


T))dt 


< 



s'{t — T)\‘^dt 





T 



0 


s''{t)fdt, 

k=l 


< 


T 

J\s'{t)\‘^dt, 
0 


(2.9) 




(s'(i) - s'(0))dt < 1 / \s''(t)?dt, 


( 2 . 10 ) 


n—1 pj' n—1 

sAt)dt 

k=0 do 


f-k + l tk 



S ioyd^dt 


< 


tk t 


n—1 


tk+1 t 


k=0 



+ (s'(0)^ ]d^dt <T [s^{t) + {s'{t)fi ]dt <{R- e)V, (2.11) 


^k ^k 


n « n—1 

/ \g\t)\^dt, \ '^rgl - / g^{t)dt < {R 
k=l i k=0 do 


efr. 


( 2 . 12 ) 


From 


■ fl2.12p it follows that 


max 


II W’l’IL?, Il[5']n|l^i) < max ( ll-sli^y^ijQ jnp 


2 

wfio 


+iR~^y^+2 j ( 2 - 13 ) 


From fl2.13p . 02.71) follows. 
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Let us know choose [u]n G V^. We simplify the notation and assume v = (s, g) = Vn{\v\n)- 
Through direct calculations we derive 


n—1 


n—1 


W|[0,T] ^ 






k=0 


k=l 


n—1 

k=0 


lt.t + + Ct, 


(2.14) 


where C is independent of r. Furthermore, we use notation C for all (possibly different) 
constants which are independent of r. By using CBS inequality we have 


^ .2 4 


-TS- 


\, = - siO)y < 




By applying Morrey inequality to s'{t) from fl2.15p it follows 

Since [n]„ G from fl2.14p . fl2.16p it follows that for all r < (2C)“^ 

IIqI|2 < (^ 

ICIIw|[o,T] — 

Therefore from fl2.14p . fl2.161) ■ fl2.17p it follows that for sufficiently small r 

n—1 n—1 n—1 

\\4wi[0,T] + + Ct 


k=0 


k=l 


k=0 


In a similar way we calculate 


n—1 


,r] + ^Slk + Ct. 


2 

W^[ 0 ,T] 


(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 


k=0 


k=l 


Hence, from fl2.18p . fl2.19p it follows that for sufficiently small r 

max (||s||^2[o,r]i II^IIw1[o,t]) < max (||[s]„|1^2, ||[^]n|l^i) + Ct, (2.20) 

From fl2.20p . fl2.8p follows. Lemma is proved. 

Corollary 2.1 (Jl^) Let either [v]n G or [v]n = Qniv) for n G Vr. Then 

\sk-Sk-i\ <C't, k = l,2,---,n (2.21) 

where C is independent ofn. 

Indeed, if n G Vr, then s' G 1F2^[0,T] and by Morrey inequality 

ll-s^llcio.T] < C'lll'S^llwjho.T] < CiR (2.22) 

and hence for the hrst component [s]„ of [n]„ = Qn{v) we have fl2.2ip . Also, if [n]„ G V^, 
then the sequence n" = Vn{[v]n) belongs to Vr+i by Lemma 1231 and the component s" of n"' 
satishes fl2.22p . Since, (s"')'(tfc) = Si^k,k = l,...,n, from fl2.22l) . fl2.2ip follows. 

Note that for the step size hi we have one of the three possibilities: hi = h, or hi = h, or 
hi < |sfc — Sfc_i| for some k. Hence, from fll.ldp and fl2.2ip . it follows that 

max hi = Oty/r), as r —)■ 0. (2.23) 

0<i<N-l 
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3 Proofs of the Main Results 


3.1 First Energy Estimate and its Consequences 

The main goal of this section to prove the following energy estimation for the discrete state 
vector. 

Theorem 3.1 For all sufficiently small r discrete state vector [M(['r']n)]n satisfies the follow¬ 
ing stability estimations: 

N—1 n N—1 


i=0 

II2 


k=l 1=0 
|2 


c[mUo,s,) + + wfwiiD) + 

n-l 

+ llx(s"(^)5 ^)IIl2(o,t) + ^ l+('Sfc+i — Sfc) ^ hiu‘^{k)^, 


(3.1) 


k=l 


where C is independent of r and 1+ be an indicator function of the positive semiaxis. 

First we prove the following lemma. 

Lemma 3.1 For all sufficiently small t, discrete state vector [M([n]„)]„ satisfies the following 
estimation: 


max 

l<k<n 


^ hiul{k) + ^ 

i=0 k=l i=0 k=l i=0 

+ l|g"lli,,„,T) + II/IIL(d) + ll7(g”(i).t)(g")'(i)llL(„.T| 

n-l 

L2{0,T) +(Sfc+i - Sfc) ^ 


(3.2) 


k=l 


where C is independent of t. 

Proof. By choosing rji = 2TUi{k) in fll.171) and by using the equality 

2TUij{k)ui{k) = uf{k) — ul{k — 1 ) + T^u^fik) 

we have 

mj — l nij — l nrij — l rrij — l 

hiuf{k) - ^ hiuj{k -1) + ffi hiulfik) + 2r hiaikU^^{k) 


i=0 


i=0 
m,' — 1 


i=0 


i=0 


2t ^ ^ hi bii^Uixi^k^Uii^k^ T Cj^Uj (fc) fik^i(.k) 


i=0 


2t [(7>”(s”)')* - xt] - 2TgtU„{k). 


(3.3) 


Using fll.bp . Cauchy inequalities with appropriately chosen e > 0, and Morrey inequality 


m-i — l 


max u^^(k) < C.||!i(i;;;)|I b.^i|„,,j < C ^ fti(M?(*:) + tiUk)) 


(3.4) 


i=0 
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where are independent of r and [M([n]„,)]„,, from fl3.3p we derive that 


m-j — 1 


m-j — 1 


mj — l 


m-j — 1 


hiUi{k) - ^ hiuj{k - 1 ) + aoT ^ hiul^{k) + hiU%{k) 


< 


i=0 


i=0 


i=0 
rrii — l 


i=0 


rrii — l 


CiT 


|(7s"(s"')0^P + |Xs"P + IS'fcP + hiU^{k) 


(3.5) 


2=0 2=0 

where Ci is independent of r. Assuming that r < Ci, from fl3.5p it follows that 

mj — l 1 mj — l 

(1 - Cir) ^ hiM-(fc) < hiM-(fc - 1) + l+(sfc - Sfc-i) ^ + 


j=0 


i=0 


Ofc_i 


CiT 


m-j — 1 


l(7..(i>")')‘P+l\7t+fe"l''+ E 


i=0 


(3.6) 


By induction we have 


mjo-l 


^ < (1 - Cir) ^ ^^(o) + ^(1 - Cir) N Cir |(7, 


j=0 


i=0 


(s")')f+ 


l=l 


ruj^-l 


1x1"^ + +l+(sz —sz_i) —1)|. 

2=0 i=mj^_^ 


(3.7) 


For arbitrary 1 < / < fc < n we have 

/ ' f ' \ —^ 

(1 - C'ir)-^+'-^ < (1 - Cir)-^ < (1 - Cit)-^ = (^ “ ^ (^.S) 

as r —)■ 0. Accordingly for sufficiently small r we have 

(1 - for 1 < / < fc < n, (3.9) 

By applying Cauchy-Bunyakovski-Schwartz (CBS) inequality from fl3.7l) - fl3.9p it follows that 

mj — l 

S /i*M?(/c) < C'2(||0 ''||L(o,so) + ll^"lli2(0,T) + Il7(s"(t),t)(s")'(t)||i2(0,r) + 

i=0 

n-1 

IL\V”(i),i)llL,„,T) + ll/llwD) + El+(*w-'’') E '“7(0). (3.10) 

Z=1 i=mj^ 
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where C 2 is independent of r. Having fl3.10p . we perform summation of fl3.5l) with respect 
to k from 1 to n and derive 

rrij^-l n mj-l n "ij-1 

^ hiul{n) + ^ hiul^{k) + hiU%{k) < 

i=0 k=l i=0 k=l i=0 

ll-#'”llL(o,«) + C3(||9"IIL,„,T) + II/IIL(D) + ll7(9”(i).t)(V)'(i)llL(0J) 

n n-1 

+ ||x(s”(t),t)||i2(0,T)+ 5^^ £ +5^1+(sfc+i-Sfc) kiuf(k) (3.11) 

k=l i=0 k=l i=mj^ 

From fl3.10p and fl3.1ip . fl3.2p follows. Lemma is proved. 

Proof of Theorem 13.11 Due to (13.21) . it is enough to show that the left hand side of (13.ip is 
bounded by the left hand side of (13.21) . By using reflective continuation property of u{x] k) 
we easily derive that 


n N—l 


^T^hiu]^{k) = 


k=l i=0 


k=l 



du{x; k) 

L 

dx 


dx < 


ps 


rsk 

du{x] k) 

L 

dx 


m-,' — 1 


dx = 2^*^TY^ hiul^{k). 


k=l i=0 


By using (I1.13P and (I2.22p we have 


(3.12) 


N-l 


N-1 


hiulik) < 2 / u^{x] k) dx h\u1^{k) < 2"*+^ f u^[x] k) dx+ 


i=0 


i=0 


Af-1 


m, — 1 


m-,' — 1 


CiT ^ hiul^{k) < 2”*+2 ^ hiul{k) + 2 ”*+^ ^ -^hlul^{k)+ 

2=0 2=0 
N-1 ”^ 9-1 N-1 

CiT ^ hiU^^lk) < 2”-+= + C^T hiufjk). 


(3.13) 


i=0 


i=0 


i=0 


From (I3.12p . (l3.13p and (13.2p . (13.ip follows. Theorem is proved. 

Let [n]„ e = 1)2,... be a sequence of discrete controls. From Lemma 1231 it follows 

that the sequence {'Pn{[v]n)} is weakly precompact in hF 2 ^[ 0 ,T] x hF 2 ^[ 0 ,T]. Assume that the 
whole sequence converges to n = {s,g) weakly in fF|[0,T] x fF2^[0,T]. This implies strong 
convergence in fF 2 ^[ 0 ,T] x L2[0,T]. Conversely, given control v = {s, g) E VJl we can choose 
a sequence of discrete controls [v\n = Qniy). Appplying Lemma l2^ twice one can easily 
establish that the sequence {Vn{\v\n} converges to n = (s,5') weakly in hF|[0,T] x 11 ^ 2 '[0,T], 
and strongly in hF2^[0, T] x L2[0, T]. In the next theorem we prove the continuous dependence 
of the family of interpolarions {u^} on this convergence. 


Theorem 3.2 Let [v]n E VJl,n = 1,2,... be a sequence of discrete controls and the sequence 
{'Pn{[v]n} converges strongly in hF 2 [ 0 ,T] x L2[0,T] to v = {s,g). Then the sequence {u'^} 
converges as r —)■ 0 weakly in IF 2 ^’ (D) to weak solution u E 1 / 2 ^’°(D) of the problem U.l]) - 
be. to the solution of the integral identity / D.idj) . Moreover, u satisfies the energy 
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estimate 


- <^(II‘?^IIl2(o,so) + llfi'llL2(0,T) + II/IIl2(d) + Il7ll^i,0(^) + llxll^i.Op)j (3.14) 

Proof. In addition to quadratic interpolation of [s]„ from fll.lip . consider two linear 
interpolations: 

s”(^) = Sfc-i + -^^(t-4-i), tk-i<t<tk,k = l,n] s'^{t) = Sn, t>T; 


Slit) = s”(t + r), 0 < t < T. 

It can be easily proved that both sequences s" and s” are equivalent to the sequence s" in 
fP2^[0,T] and converge to s strongly in 11^2^ [0,T]. In particular, 


sup||Si||h^i[o,t] 


<a 


(3.15) 


where (P* is independent of n. 

Our next goal is to absorb the last term on the right hand side of fl3.1l) into the left hand 
side. We have 


71—1 


"^4+1-1 


^l+(sfc+i-Sfc) hiulik)< 


k=l 


2X]l+(sfc+i - Sfe) / u^ix;k)dx+-Y'^+i^k+i-Sk) Y (3-16) 

k=l k=l i=mi 


Note that if Sk+i > Sk, then all the factors hi in the second term are bounded by — Sk 
and by using fl2.2ip we have 


71—1 




71—1 


Y i+(sfc+i - Sk) Y ^ ^ / 

i=mi k=l 


Sk + l 


k=l 


du 


dx 


dx 


(3.17) 


Due to reflective continuation we have 


fSk + l 

du 

^dx < 2”*“^ 

j-Sk 

du 

' Sk 

dx 


Jo 

dx 


771,-1 


dx = 2"‘ ^Y hi^ixik)- 


i=0 


From (I3.17P and (I3.18p it follows that 

71—1 




n-1 ruj-l 

Y l+(sfe+i - Sk) Y < 2”*-^(C")^^ ^ £ hiul^ik) 

k=l i=mi k=l i=0 


(3.18) 


(3.19) 


Assuming that r is sufficiently small and by using fl3.16p - fl3.19p in fl3.ip . we absorb the last 
term on the right hand side of fl3.19l) into the left hand side of fl3.1l) and derive modified 
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fl3.ip with a new constant C\ 


N-l 


n N—1 


“ “ i=0 k=l i=0 

C {\\^'^\\l2{0,so) + llfi'"'llL(0,T) + ll/llL(0) + ll7('S”(t),t)(s"')'(t)||^2(0^r) 

r^k+i s 

+ llx(s”(t),t)||i 2 (o,T) + 5^1+(sfc+i - Sfc) / u‘^{x;k)dxV 

k=l 

We can now estimate the last term on the right hand side of fl3.20p as in [1]: 

/•-'ii-1 1 n—1 fifc+l 


(3.20) 


^ rsk+i " ^ rH 

^l+(sfc+i-Sfc) / U^{x-,k)dx = '^l+{sk+i-Sk) / 

k=l k=l 


s”) (t)u^it)] k)dt = 


fc=l 
^fc+1 


^l+(sfc+i-Sfc) f - r)) dt = 

k=i ^ ' 

r^k , s 2 

5^ l+(Si+i - Si) / (i^) w{y{i^(t),t)j dt. 


(3.21) 


By applying CBS inequality we have 


^ ( r^k+i 

^l+(sfc+i -Sfc) / U^{x]k)dx < ||(s)^)'|U2[0,T]||M^(si(^),^)||i4[0,T]- (3-22) 

k=l 


From the results on traces of the elements of space V 2 {D) ([211 El EH]) it follows that for 
arbitrary u G V 2 {D) the following inequality is valid 

||-u(Si(f),f)||L4[o,T] <C\\u\\v2{D)i (3.23) 


with the constant C being independent of u as well as n. From (13.151) . (13.221) and (13.231) it 
follows that 


^ l+(sfc+i - Sfc) / U^{x]k)dx 

k=l ^k 

If the constant C* from fl3.15l) satisfies the condition 


< C.C\\u^ 


2 

V2{D)- 


(3.24) 


a < {CC)-^ 


(3.25) 


then from fl3.20l) and fl3.24p it follows that 

S C'(ll'#’”llis(0,«) + lls”lli 2 ( 0 ,T) + II/IIlj(D) + 

||7(s”(«).()(s”)'(«)||i4„^„ + ||x(s”(«),«)ffa(„,„). (3.26) 


where C is another constant independent of n. By applying the results on the traces of 
elements of Wl’^{D) ([H EHj) on smooth curve x = s^(t), Morrey inequality for (s"')' and 
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we have 


||7(s"'(t),t)(s"')'(t)||2,2(o,r) < II(s"')1Ic[o,t]||7(s"'(^);^)||l2[o,t] < ^7,\\i\\wI’°{d) 

\\x{s^{t)ML,[o,T]<C,\\x\\w..o^r>)^ (3.27) 

where C 3 is independent of 7,7 and n. Hence, from 03.261) . 03.27P it follows the estimation 




c 110 " 


|2 

1^2(0,so) 


+ 11 ^/^ 


|2 

\L2{0,T) 


+ 


2 

L 2 (D) 



llxll 




(3.28) 


with C being independent of n. 

If 03.25P is not satished, then dne to 02.2ip we can partition [0,T] into hnitely many 
segments [tnj_i,tnj], 7 = 1,? with tno = 0, tn^ = T in snch a way that by replacing [0,T] 
with any of the snbsegments [tnj_i,tnj] fl3.15p will be satished with C* small enough to obey 
03.25p . Hence, we divide D into hnitely many subsets 


= Dn < t < t^.} 

such that every norm ||M"||y 2 (n») uniformly bounded through the right-hand side of 03.28p . 
Summation with j = 1,..., g implies 03.28p . 

Since (0"',fi'”') converge to ( 0 ,( 7 ) strongly in L2[0,so] x L2[0,T], from 03.28P it follows 
that the sequence {m"|- is weakly precompact in Let u G be a weak 

limit point of u" in H 2 ’°(Z 7 ), and assume that whole sequence {«"} converges to u weakly 
in Wl’^{D). Let us prove that in fact u satishes the integral identity fll.lOp for arbitrary 
test function $ G H^ 2 ^’^(H) such that ~ density of C^(f2) in H^ 2 ^’^(f 2 ) 

it is enough to assume d* G C^(H). Without loss of generality we can also assume that 
<h G C^{Dt+t), *h = 0, for T < t < T -I- r, where 


Dt+t = : 0 < X < I + 1, 0 < t <T + t} 


Otherwise, we can continue <h to Dt+t with the described properties. Let 
^i{k) = <h(a;i, 4), A: = 0, • • • , n -b 1, i = 0, • • • , iV 


and 

<h"(a;,t) = <hi(/c), = ^ix{k), = ^it{k + 1 ), for 4 _i <t<tk,Xi<x < Xi+i. 

Obviously, the sequences {$"}, and {$[} converge as r —)■ 0 uniformly in Z) to 
and ^ respectively. By choosing in fll.l7p 7 * = T^i{k), after summation with respect to 
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k = l,n and transformation of the time difference term as follows 


rrin—l 


n—1 ^Jk+i ^ 


hiUii{k)^i{k) = -y^ T hiUi{k)^ii{k + l) - y^ hiUi{ 0 )^i{l) + 


k=l i=0 
n—1 


k=l i=0 
Pk-l 


i=0 


" ^ rtk rsk+1 

yy sign{sk - Sk+i) hiUi{k)^i{k) = - 

k=l i=ak k=l Jo 


dx dt- 


-•si 


n—1 


0 k-i 


0”(a:)$^(x, T)dx + '^ sign{sk - Sk+i) 

k=l 'i'=0'k 

l-T ^s{t) nsi 

' 0"(x)$'^(x, r) dx - 




'Xi 


f*T-7 


dx dt 


>0 Jo 


" ^ rtk rsk+i ^ ^ ^ ^ 

yy / / dxdf - - '^sign{sk - s^+i) 

k=l J^k-l Js{t) J 


m(x; /c) — Uix{k){x — Xj) j$j(A:)dx = 

(s7)'(f)n-((.7)(f),f)$-((s7)(f),f)df 
Pk-l 

(3.29) 


where 


Offc = ), ( 3 k = max(mj,^,mj ), 


we derive that 


'0 JO 


r-so 


T rs{t) r fly'^ Su^ 

a—^l - &—}> r) dx 

ox ox ' ' 


r-r-7 


s7)'(t)n"((5")(t),f)«h-((s7)(t),f)df+ / <7"(t)«h-(0,f)df 


rT ^ 


+ 


7(s"(f), f)(s")'(f) - x(s”(t), t)) <h-(s"(t), t) dt-R = 0 


(3.30) 


where 

Ok fSk 

«=E/ / 

k = l ^'^k -1 ’J s(t) 
C^k r^k 


du 


du 


+-f*' 1 * * - E 


dx 


n—1 




dx dt 


't'^k — l ^ s(t) 


k=l J^k-l Js-^it) 

Ok rsk +1 




7(s”(f),f)(s'^)'(f)-x(s''(t),t))J-^dxdf+ / (/)”(x)<h^(x,r)dx 


'so 


^ ( Ok j-Sk + l ^ ^ 

y^, / / M^<h[dxdt - - y^sz5fn(sfc - Sfc+i) dyia;(A:)<hi(A;) 

fc=l “'ife-i ^ fc=l i=aj. 

First note that the sequence {n"^} is equivalent to the sequence {u^} in strong, and accord¬ 
ingly also in a weak topology of L 2 {D), and hence converges to u weakly in L 2 {D). Indeed, 
by using fl3.ip we have 


n N—1 


1“^ “^1112(15) 3 


- r hiul^{k) max hj 0, 


as n —)■ cx). 


(3.31) 


k=l i=0 
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Let 

n 

^ = U : 4-1 <t <tk, mm{s{t),Sk) <x < max(s(t), s^)} 

k=l 

|A| denotes the Lebesgue measure of A. Since s"‘{tk) = Sk, we have 


|A| < ^ / / \s'{T)\dTdt + '^T\s{tk) - s"-{tk)\ < 

k=l k=l 

V^||s 1 |l 2 ( 0 ,t)T + T'IIs - s’^||c[ 0 ,r] -t 0 as r -)■ 0 

and all of the integrands are uniformly bounded in L^{D), it follows that the hrst term in 
the expression of R converges to zero as r —)■ 0. In a similar way one can see that the second, 
third and hfth terms also converge to zero as r —)■ 0. The fourth term in the expression of 
R converges to zero due to Corollary 12.11 and uniform convergence of {*h^} in D. To prove 
the convergence to zero of the last term of R, let 


n—1 

^ = U <tk,dk= min(sfc, s^+i) < a; < 4+i = max(sfc, Sfc+i)} 

k=l 

From Corollary 12.11 it follows that 


Since 


we have 


A| < Ct —>-0, as r —)■ 0. 


/3fc-l 

E*- 

i=0‘k 


^k '^fc+1 


n—1 


/3fc-l 


n-1 


Y,sign{sk - Sk+i) Y, hlui,{k)^i{k) < Y I 

du’ 


f'ife+i 


k=l 


i=ak 
rdk+i 

T / 


k=l 


dx 


\^'^\dx < 


k=l 

du’ 




du’ 


dx 


\^'^\dx < 


dx 


L2{A) 


Id)" 


Il2(a) 


(3.32) 


Since the integrands are uniformly bounded in L 2 {D), the expression in fl3.32p converges to 
zero as r —)• 0. Hence, we have 

lim i? = 0 (3.33) 

r—>-0 

Due to weak convergence of to u in weak convergence of iR to u in L 2 {D) and 

uniform convergence of the sequences {$""}, and {d)[} to d), || and ^ respectively, 

passing to limit as r —)■ 0, it follows that hrst, second and fourth integrals on the left-hand 
side of fl3.30p converge to similar integrals with n"" (or u^), d>^, d)[, d)'"(x, r), g^(t), cfR^x) and 
d>'"(0,f) replaced by M,d>, d)(a;,0), g{t), 4>{x) and d>(0,f) respectively. Since s" converges 

to s strongly in 1 T 2 ^[ 0 ,T], the traces 7 (s""(f), (t)), x(s”(t),f) converge strongly in L 2 [ 0 ,T] to 
traces 7 (s(f), (f)), x(s(t),f) respectively. Since d)'"(s”(t), t) converge uniformly on [0,T] to 
d)(s(t),t), passing to the limit as r —)■ 0, the last integral on the left-hand side of fl3.3Up 
converge to similar integral with s" and d*"" replaced by s and d*. 
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It only remains to prove that 


pT—r pT 

lim / t) hi = / s'(t)u(s(t),t)^(s(t),t) (3.34) 

Jo 

Since {s”} converges to s strongly in hh 2 ^[ 0 ,T], from fl3.24p it follows that is 

nniformly bounded in L 2 [ 0 ,T] and 

||M^(s"(t),t) - M’‘(s(t),f)||L 2 [ 0 ,T]-t 0 as r0 (3.35) 

Since {«’’} converges to u weakly in it follows that 

u'^{s{t),t) —)■ u{s{t),t), weakly in L2[0,T] (3.36) 

Since {<I)^(s”(f), f)} converges to uniformly in [0, T], from fl3.35p . fl3.36p . fl3.34p easily 

follows. 

Passing to the limit as r —)■ 0, from fl3.30p it follows that u satishes integral identity fll.lOp . 
i.e. it is a weak solution of the problem fll.ip - fll.4p . Since this solution is unique (|2l]) h 
follows that indeed the whole sequence {u^} converges to n G f 2 ^’°(fl) weakly in fP 2 ^’°(fI). 
From the property of weak convergence and fl3.27|) . fl3.14|) follows. Theorem is proved. 

In particular, Theorem 13.21 implies the following well-known existence result i |24] ): 

Corollary 3.1 For arbitrary v = {s,g) G Vr there exists a weak solution u G ^^’^(fl) of the 
problem U.l\) - ^1.4\ ) which satisfy the energy estimate ( [J. j^[ ) 

Remark; All the proofs in this section can be pursued by using weaker assumptions 

0 G L2[0, /], 7 , X e W^’\D), a G L^{D), 

and fll.6p instead of conditions imposed in Section 11.51 The only difference would be to 
dehne (fi as a Steklov average 

r^i+\- 

(j)i = -r (j){x)dx, i = 0,...,N — l 
hi Jxi 

and replace the norm of 0"" in the hrst energy estimate through norm of 0. 

3.2 Second Energy Estimate and its Consequences 

Let given discrete control vector [n]„, along with discrete state vector [n([n]n)]„, the vector 

[«(Hn)]n = (m(0), m(1), ..., u{n)) 


is dehned as 


Ui{k) 


Ui{k) 0 < i 

Urrij {k) mj < i < N ,k = n. 


The main goal of this section to prove the following energy estimation for the vector n([n]ri,). 
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Theorem 3.3 For all sufficiently small r discrete state vector [^(Hn)]^ satisfies the follow¬ 
ing stability estimation: 


n I'tj 


mn — l 


n I'f'j 


rriA —1 




< 


c 


110" 


i=0 k=l i=0 

2 , II ni|2 


|2 

l-J^2[0,So] 


+ ll0llvyi[o,so] + 11^ 


'wfio,T] 


k=l i=0 




Wf [0,T] 


+ ||x(s"(t),t)||' 1 +"""' 

W}[0,T] 


L2{D) 


(3.37) 


Proof: Note that if Sk-i > Sk then ufik) can be replaced through ufik) in all terms of 
fll.lTp . By choosing rji = 2TUii{k) in fll.lTp and by using the equality 

‘^raikUix{k)Ui^t{k) = OikU^ffk) - ai^k-iujx{k - 1) - raiktU^xi^ - 1) + T^aikU^ikti^)^ {3.3S) 


we have 


rrij — l rrij — l rrij — l rrij — l 

hiUikulxik) - 'Yl hiOkk-iulxik - 1) + 2r Kulfik) + Y 
2=0 2 = 0 2 = 0 2=0 

VTij — l rrij — l rrij — l 

= r Y hiaiktulx{k - 1) + 2r hihikUix{k)Uit{k) + 2 tY^ hidkufikffffk) 

2=0 2 = 0 2 = 0 

nij-l 

-2t Y hifikUffik) - 2t “ XsA - ‘^rg^Uoffk) (3.39) 

i=0 


If Sk -1 < Sk, then ufik) can be replaced through Ui{k) in all but in the term including 
backward discrete time derivative in fll.171) . The latter will be estimated with the help of 
the following inequality: 

TTij — l rrij — l rrij — l 

2t Y hiUiffk)uffk) >'rY Y “ 1 ) ( 3 . 40 ) 

2=0 2=0 i=mj, ^ 


To prove fl3.40p . we transform the left hand side with the help of the CBS and Cauchy 
inequalty with e = r to derive 


mj — l 


2—1 


2r Y hiUu{k)uft{k) = 2 t Y Y Y “ 1) 

p=^ik-i 


2 = 0 


2 = 0 


rrin — l 


mn — l 


'"dk-l 
2 — 1 


>rY hiUffik) - - Y Y “ 1) 


2=0 
m-j — l 


>rY hiul{k) - 


mj — l 

|2 


^k —1 


Y1 


(3.41) 


2 = 0 


2=m 


‘dk-i 
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which implies fl3.40p due to fl2.2ip . Hence, in general fl3.39p is replaced with the inequality 


777.0—1 


m-j—1 


rrin—l 


777.0—1 


^ ^ hitti^k-iui^ik - 1) + r ^ hiU%{k) + hiaikuf^j{k) 


i=0 


i=0 


i=0 


i=0 


rrii — l 


m, — 1 


m, — 1 


< 


r ^ hittiktuf^ik - 1) + 2r ^ hibikUi^{k)Uii{k) + 2r ^ hiCikUi{k)uft{k) 


1=0 


i=0 


i=0 


rrii — l 


m, — 1 


+ {C'fr ^ hiul^{k - 1) - 2r ^ hifikUft{k) 




7=0 


-2r - Xs"] - 2 rc/>o 7 (^) (3-42) 

By adding inequalities fl3.42p with respect to k from 1 to arbitrary p < n we derive 

rrijp-l p rrij-l p rrij-l 

h,a,puUp) + TY, hiuljik) + hittikU-^^^k) 


i=0 


k=l i=0 
rrii — l 


k=l i=0 

p mj-1 


< {C'f^ ^ l + (Sfc - Sk-l) ^ ^i^iktulAk - 1) 

k=l k=l 7=0 

p rrij — l rrij — l rrij — l 

+2r E[E hibikUi^{k)ua{k) + 2 ^ hiCikUi{k)u{i{k) - 2 ^ hJikUaik) 

k=l i=0 i=0 i=0 

mj^—l p p 

+ - 2r [(7,n(s”)')^ - j(/c) - 2r^^^no,7(^) 


(3.43) 


i=0 


k=l 


k=l 


By using fll.6l) and by applying Cauchy inequalities with appropriately chosen e > 0, from 
fl3.43p it follows that 


mjp-l 


Oo 


i=0 


p ruj 


rriA —1 


P ^3 


777o — 1 


P ^3 


rrii — l 


h^uUp) + 7^J2Yl hiU%{k) + OoT^ EE hiul^t{k) < T 


k=l i=0 

n ^rrii-l 


k=l i=0 


rrii — l 


+Ct EE hiuf^k) + ^ ^ hiuf^[k) + E'^./ii +C’ 5: hi(j): 


k=l i=0 
n 


i=0 


i=0 


k=l i=0 




i=0 


+2r^|(7sn(s")')^-xJn| \um,:t{k) \ +2r^|^^| 1 mo,7(^)| (3-44) 

k=l k=l 

where C is independent of n. First term on the right hand side will be estimated as follows: 

p rrij — l n m, —1 


r 


^ _■> "■ / 1 rxi+x ptk r 

^ hittiktuUk = / / 

k=l i=0 k=l i=0 '^7-1 Jt- 

da{x, t) 


1 /■*'= /■* da{x,i) 

~W~ 


d^dtdxu^^{k — 1 ) 


< 2 / esssupo<x<i 
Jo 


dt 


777o — 1 


hiul{k) + C hi(t). 


(3.45) 


i=0 


i=0 
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Due to arbitraricity of p, from fl3.44p it follows 


n ruj 


rrin — l 


n 


m-i — l 


+ 2 YY Y Y 

i=0 k=l i=0 k=l i=0 


<2 esssupo<^<i 


k=l i=0 

da{x, t) 


dt 


m-,' — 1 




dt max > hiujjk) + (7 > hx 

l<k<n ^ ^ 


I 


i=0 


i=0 


_ rrii — l 


+^'^Y Y + Y + Y 

k = l f 2 = 0 2=0 2 = 0 

n n 

+2ry] |(7,n(s”)y - X^n| \urn^;t{k) \ +2r^|c/^| |Mo,t(^)| 


k=l 


k=l 


If 


CSSSUpQ^x^l 


da{x, t) 


dt 


dt <C Oq 


(3.46) 


(3.47) 


then the hrst term on the right hand side of fl3.46p is absorbed into the hrst term on the 
left hand side. If fl3.47p is not satished, then we can partition [0,T] into hnitely many 
subsegments which obey fl3.47p . absorb hrst term on the right hand side into the left hand 
side in each subsegment and through summation achieve the same for (I3.46p in general. 
Hence we have 


m-j — 1 


n ’<^3 


m-i — l 


n >'^3 


m, — 1 


max 

l<fc<r2 


,2 

IX 


(3.48) 


2=0 k=l 2=0 k=l 2=0 

^rrij — l TTij — l rrij — l rrij^—l 

<ctY, E '><“?('=) + E '>■“?.('=) + E +c E 

k=l f 2 = 0 2=0 2=0 2=0 

n n 

+CtE |(7.»(s”)r - xtl + CrE Wl |fio,!(*;)| 

k=l k=l 

with some C independent of n. 

Since 7 , x G WP{D) we have 7 (s’^(f), t), x(s”(f), t) G 1 H 2 ‘‘[ 0 ,T] (|28l El |2l]) and 
||7(s (t),2)||^^i^^^j < C||7||^yi,ip), ||x(s (2))^)X|jp^j < C'llxllwj.hD)’ 

where C is independent of n. According to Lemma 1231 V^{[v\^'] G Vr+i. By applying Morrey 

1 

inequality to (s”)' we easily deduce that 7 (s"'(f),f)(s"')'(f) G 1H2'‘[0,T] and moreover, 


(3.49) 


h{s-{tit){s-y 


w} [o,r] 


< C'i||7(s"'(2),t)||^i^^^j||s"||w|[o,T] < C'||7||^i,i(^), (3.50) 


where C is independent of n. 

Let w{x,t) be a function in W^’^iD) such that 

w{x,0) = (j){x) for a; G [0, So], a{0,t)wx{0,t) = g'^{t), for a.e. f G [0, T] 
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(3.51) 
























a{s'^{t),t)w^{s'^{t),t) = 7 (s’^(t),t)(s”)'(t) - x(s”(t),t) for a.e. t e [0,T] 


(3.52) 


and 






< C 






+ ||7(.-(t),t)(.'^)'(t)-x(^"(t),t)ll 

W 2 [0,-/ J 


(3.53) 


The existence of w follows from the resnlt on traces of Sobolev fnnctions [5l |28]. For example, 
w can be constrncted as a solntion from of the heat equation in 


= {0 < a; < s'^(t),0 < t < T} 

under initial-boundary conditions fl3.5ip . fl3.52p with subsequent continuation to W 2 '^{D) with 
norm preservation [33l IM] . 

Hence, by replacing in the original problem fll.ip - fll.4p u with u—w we can derive modihed 
fl3.48p without the last three terms on the right-hand side and with /, replaced by 

F = f + wt- {aw^)^ - bw^ -cw e L 2 {D). (3.54) 


By using the stability estimation fl3.2p . from modihed fl3.48p . fl3.53p and fl3.54p . the following 
estimation follows: 


777.0—1 


n 


777,0—1 


n fflj 


rriA —1 


max 

l<k<n 


hiuUk)+^YY 


< 


i=0 


k=l i=0 


k=l i=0 


c 




I n2[o,so] 






'W2^[0,T] 


71—1 


”^4+1-1 




2 

L2{D) 


^l+(sfc+i-Sfc) Y hiul{k) 


k=l 


(3.55) 


By estimating the last term on the right hand side of fl3.55p as in the proof of Theorem 13.21 
fl3.37p follows. Theorem is proved. 

Second energy estimate fl3.37p allows to strengthen the result of Theorem 13.21 

Theorem 3.4 Let [n]„ G V^, n = 1, 2,... he a sequence of discrete controls and the sequence 
{Vn{[v]n} converges weakly in hF 2 ^[ 0 ,T] x hF 2 ^[ 0 ,T] to v = {s,g) (i.e. strongly in hF 2 ^[ 0 ,T] x 
L 2 [ 0 ,T]) to V = {s,g) for any 5 > 0, 


n {x < s{t) — 6 , 0 < t < T} 

Then the sequence {h'^} converges as r —)■ 0 weakly in Wl’^iVt') to weak solution u G Wl’^iVt) 
of the problem i.e. to the solution of the integral identity M.tA) . Moreover, u 

satisfies the energy estimate 


\u\ 


wf\n) 


<c{U\ 


|2 

lwi(0,so) 


+ 


2 

1 

wL 


[0,T] 


+ 


2 

L2{D) 


+ 


wf\D) 


+ 


wf\D) 


(3.56) 


Proof; Let 1 0 be arbitrary sequence and 


= {(a^, t) : 0 < X < s{t) — Cm, 0 < t < T} 
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Note that the sequence s" converges to s uniformly in [0,T]. Calculate 


I 

1“ 


E 


ftk /■s(t)-e„ 


+ 


k = l -^0 


dv7 


dx 


+ 


du" 


dt 


dx dt 


Denote s™ = a:-, where 


i = max i : max s{t) — < Xi < max s{t) - 


and substitute to derive 


rtk 

k=iJtk-iJo 


\u{x; k — 1) + Ui{x; k){t — tfc-i) dx dt+ 


r-4 


E 

k = l Jo 


du{x; k — 1) dui{x; k) 


dx 


+ 


dx 


-{t-tk-i) 


dx dt+ 


Hk r^l 


E 

k=l J^k-l Jo 


\ui{x-, k)\‘^ dxdt 


From which by using fll.lSp it follows that 

n i—1 


n z—1 


I -r||2 
1“ 


< c^T Y ^ “ i)+ 


k=l 0 
n i—1 


k=l 0 
n i—1 


+ tYJ 2 + 


'^ixd 


k=l 0 


k=l 0 


(3.57) 


where C is independent of n or r. Our goal is to prove that the right hand side of fl3.57l) is 
bounded by the left hand side of fl3.37p for sufficiently large n. It is sufficient to prove the 
following claim: for fixed e^, there exists N = N{em) such that for\/n > N 


sff < min(sfc, Sk-i), k = l,...,n 
Indeed from (I3.58p it follows that 

n i—1 n i—1 n i—1 


(3.58) 


T 

k=l 0 
n mj-l 


k=l 0 
n riij-l 


k=l 0 
n mj-l 


< ^ 5 ^ hiul^{k - 1 ) + ^ hiulik) + r'^YYl 

k=l 0 k=l 0 k=l 0 

To prove fl3.58|) . we hrst show that for sufficiently large n and all tfc-i <t<tk 


sit) -Sk< — 


(3.59) 


(3.60) 


We have 


sit) - Sk = sit) - s(4) + s(4) - s^itk) + s"(4) - Sk 
E ||s1|c[o,r]'?‘ + ||s — s'^||c[o,r] + s'^itk) — Sk 
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Observe that 


so 


s”(4) 


Sk + Sfc-1 
2 


s{t) — Sfcl < 11 s'11 (7(0,T] '7" + ||s — s”||c[0,T] + 


Sfc—1 S/j 

2 


By Money’s inequality and fl2.ip . 


|s(t) — Sk\ < ^OIIs'lliy^qo^T] + 7- + 


lc[o,r] 


(3.61) 


By uniform convergence of s"" —)■ s on [0, T] there exists Ni = iVi(em) for which fl3.60p holds. 
Similarly, there exists some N 2 = N 2 {em), such that for n > N 2 and for all tk-i < t < tk, 
fl3.60p is true with Sk replaced by Sk-i- Hence, fl3.58p follows with N = max(A^i, N 2 ). 

Applying fl3.57p . fl3.59p . second energy estimate fl3.37p . and the hrst energy estimate fl3.ip . 
derive 




110" 


|2 

Il2[0,So] 


+ II0II 


Wi[0,so] 


+ 


L2{D) 


+ 




n—1 


'•Jfc+1 


-1 


+11x11 


W^’\d) 


..77 II 2 


+ 111? \\wl^\o,T] 


^^l+(sfc+i Sfc) 'y ^ hiUj^{k) 


k=l 


(3.62) 


By estimating the last term on the right hand side of fl3.62p as in the proof of Theorem 13.21 
we derive 


\u 


||2 



(3.63) 


Since 0” — )■ 0 strongly in + 2 ( 0 , sq], g'^ ^ g weakly in H0^[O, T], the right hand side is uniformly 
bounded independent of n. Hence, is weakly precompact in It follows that 

it is strongly precompact in L 2 {klm)- Let m be a weak limit point of {if} in and 

therefore a strong limit point in L 2 (Hm)- From anothe side the sequences {u^} and {u^} are 
equivalent in strong topology of L 2 (Hm)- Indeed, we have for all n > N{rn) 


\u 


u 


|2 

I r.2 (O7] 


) < 2r- 


n rrij-l 


2 ^ 2^ \hiUit{k) + =0(r), as r ^ 0, (3.64) 


k=l i=0 


due to second energy estimate (I3.37p . Therefore, m is a strong limit point of the sequence 
{u^} in L 2 (Hm)- By Theorem 13.21 whole sequence {u'^} converges weakly in hF2^’°(r2) to the 
unique weak solution from 1 / 2 ^’°(H) of the problem fll.ip - fll.4l) . Hence, m is a weak solution 
of the problem fll.ip - fll.4p and we conclude that whole sequence {if} converges weakly in 
W 2 ^{Vtra) to M G {0.^) which is a weak solution of the problem fll.ip - fll.4p from 14^’°(r2). 
Hence, Ut exists in kirn and ||Mt||L 2 (o^) is uniformly bounded by the right hand side of fl3.63p . 
It easily follows that the weak derivative Ut exists in H, and u G W^'^ikl). By using well- 
known property of the weak convergence, and passing to limit hrst as n -foo, and then 
as m —)■ - 1 - 00 , from (I3.63p . (13.561) follows. Theorem is proved. 

In particular. Theorem 13.41 implies the following existence result: 
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Corollary 3.2 For arbitrary v = {s, g) ^ Vr there exists a weak solution u G of the 

problem U.l\) - ^1.4\ ) which satisfy the energy estimate HU.56]) . By Sobolev extension theorem 
u can be continued to with the norm preservation: 


\u\ 


wf\D) 


<c m 


\2 

lvyi(o,so) 


+ 




[0,T] 


+ 


LiiP) 


+ 


wf\D) 


+ 


wf\D) 


(3.65) 


Remark: In fact, we proved slightly higher regularity of m, and both in Theorem 13.41 and 
Corollary 13.21 or on the left-hand sides of fl3.56p or 03.651) can be 

replaced with 


\u\ 


= max 

0<t<T 


\U 


^)llw2ho,s(t)] 


+ \\ut 


|2 

\L2{n) 


or M 


= max 

0<t<T 




The proof of the Theorem II. II coincides with the proof of identical Theorem 1 1.1 1 in [T]. The 
main idea is that hrst and second energy estimates imply weak continuity of the functional 
x7{v) in 11^2^ [0,T] X Wff)^T]. Since Vr is weakly compact existence of the optimal control 
follows from Weierstrass theorem in weak topology. 

We split the remainder of the proof of Theorem 11.21 into three lemmas. 

Lemma 3.2 Let jr*(±e) = inf fjfv), e> Q. Then 

W±e 


hmj7*(e) = J^ = hmX(-e) (3.66) 

The proof of Lemma [3.21 coincides with the proof of identical lemma 3.9 from [T]. 

Lemma 3.3 For arbitrary v = {s,g) G Vr, 

lim In{Qniv)) = J{v) (3.67) 

n—>-OD 


Proof: Let v G Vr, u = u{x,t;v), Qn{v) = [n]„ and [M([w]n)]n be a corresponding discrete 
state vector. In Theorem 13.41 it is proved that the sequence {vT} converges to u weakly 
in for any hxed m. This implies that the sequences of traces {'u’'(0,f)} and 

{u'^{s{t) — em,t)} converge strongly in L 2 [ 0 ,T] to corresponding traces u{0,t) and u{s{t) — 
Cm, t). Let us prove that that the sequences of traces {m"^( 0, f)} and {u'^{s(t)—em, t)} converge 
strongly in L^[0,T] to traces u{0,t) and u{s(t) — em,t) respectively. By Sobolev embedding 
theorem ([5l [28]) it is enough to prove that the sequences {u'^} and are equivalent in 
strong topology of Wl’^iVLrn)- In Theorem 13.41 it is proved that they are equivalent in strong 
topology of L 2 {Llra). It remains only to demonstrate that the sequences of derivatives ^ 
and are equivalent in strong topology of L 2 (flm)- Following the proof of the Theorem 13.41 
from the second energy estimate fl3.37p it follows that for ah n > N{m) 


du'^ 

diF 

dx 

dx 


2 

L 2 (^m ) 



n ^ 3 -^ 




as r —)■ 0. 


(3.68) 


Let V{t) = tV{t) = pf ,ii tk-i < t < tk, k = 1,... ,n. We have 


~ ^||l2[o,t] 0, Wp'^ — p\\l^[o,t] —)■ 0 as r —)■ 0 


(3.69) 
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We estimate the first term in X„(Q„(r;)) as follows 

n n 


/Sqt ^ |Mo(fc) — =/So ^ / \uo{k) — dt =/3o / |m'^( 0, t) — W(t)p (it (3.70) 

j=i k=i Sifc-i So 


From fl3.69p it follows that 

n 

lim P„T ^ |ti(,(4) - = A||ti(0, t) - 


(3.71) 


k=l 


We estimate the second term in X„(Q„(n)) as follows 


k=l 


+AX]/ |M(g(t); k) - dt + /3i 


A:=l -tife-l 

tfc / r^k Qyi 


—/i’'(t)) dxdt 


k=l 


— ^^k — 1 s(t) 


dx 


dx I dt — Ii 12 -\-13 (3.72) 


Since 


(9ii^ 


dx 


L 2 (D) 


and ||M'^(s(t), t) — /n'^||L 2 [o,T] are uniformly bounded, and {s”} converges to 


s uniformly on [0,T], by applying CBS inequality and fl3.6ip it easily follows that 

lim Ji = 0, lim I 3 = 0 

n—^oo n—>-oo 

It remains to prove that 


(3.73) 


lim /2 = lim/di / \u'^{s{t),t) — jx'^{t)f dt =/3i / \u{s{t),t) — ix{t)f dt (3.74) 


n^oo n—^00 


Since /i"^ —)■ /r strongly in in X 2 [ 0 ,T] it is enough to show that 

\\u'^{s{t),t) - M(s(t),t)||i 2 [o,T] -t 0 
as r ^ 0. For any hxed m > 0, estimate 


(3.75) 


\\u'^{s{t),t) -M(s(t),t)||L2[0,r] < \\u'^{s{t),t) -u^{s{t) - em,t)\\L^[0^T] + 

+ ||M^(s(t) - em,t) - M^(s(t) - em,t)|U2[0,T] + 

+ ||M^(s(t) -e^,t) -u{s(t) - em,t)||L2[0,T] + \\u{s{t) - e^,t) - M(s(t), t)||L2[o,r] (3.76) 


Estimate the hrst term on the right-hand side of fl3.76l) as 

\\u'^{s{t),t) - M^(s(i) - e,n,i)IU 2 [ 0 ,T] = 

By CBS inequality and the hrst energy estimate 


du^{x,t) 2 


' s{t)-€„ 


dx 


dx 


\ 1/2 


dtj 


(3.77) 


|M^(s(t),t) - M^(s(t) - em,t)||L2[o,T] < \/e 


du'' 


dx 


L2{D) 


< c.re 


(3.78) 
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for C independent of n. Similarly, the last term in fl3.76p is estimated by using CBS and 
energy estimate (13.1411 : 


\u{s{t) - e^,t) - u{s{t),t)\\L^io,T] < 


(3.79) 


Fix e > 0 and find M such that for all m > M, C^/e^ < e/4. Taking m = M, it follows 
from dSTnD-dSZi 

||M^(s(t),t) - u{s{t),t)\\L2[0,T] <1 + ll“^(■5(^) -eM,t) -M^(s(t) - Cm, t) ||l2[0,T] 

+ ||h^(s(t) - eM,t) - u{s{t) - eM,t)||L 2 [ 0 ,T] (3.80) 

The second term in (13.8011 is estimated through Sobolev embedding of traces as 

||n^(s(t) - Cm, t) - u^{s{t) - Cm, t)||L 2 [o,r] < C\\u^ - (3.81) 

By (13.6411 . (I3.68P there exists Tq{M) > 0 such that Vr < tq 

||M^(s(t) -eM,t) -M^(s(t) - eM,t)\\L2[o,T] < ^ (3.82) 

According to well-known compact embedding theorem weak convergence of ■u’’ —)■ u in 
implies strong convergence of traces ir\x=s{t)-eM '^\x=s(t)-eM -^ 2 ( 0 ,T] 

that is, there exists ti(M) such that for all r < ti we have 


\u'^{s{t) - eM,t) - U{s{t) - eM,t)\\L 2 [ 0 ,T] < I 


(3.83) 


Hence, according to (13.8011 . (13.8211 . (I3.83p . for arbitrary e > 0 we can hnd T 2 = min(ro;ri) 
such that for all r < r 2 we have 


\\u^{s{t),t) -M(s(t),t)||L 2 [ 0 ,T] < e, 
which proves (13.7411 . Lemma is proved. 

Lemma 3.4 For arbitrary [n]„ G 


(3.84) 


lim (j{Vn{\v]n)) -Tn{\v] 


= 0 


(3.85) 


Proof: Let [v]n G and n” = = Vni[v]n)- From Lemma [2^ it follows that the 

sequence {Vn{\v\n} is weakly precompact in fF 2 ^[ 0 ,T] x hF 2 ^[ 0 ,T]. Assume that the whole 
sequence converges to h = {s,g) weakly in 1F2^[0,T] x hF2^[0,T]. This implies the strong 
convegence in 1F2^[0,T] x L2[0,T]. From the well-known property of weak convergence it 
follows that V G Vr. In particular s” converges to s uniformly on [0,T] and we have 


lim max \s"'(ti) — s(tj)| = 0 
n—^-oo 0<2<n 

We have 

MM - = MM - 

Since 77’(n) is weakly continuous in 1F2^[0,T] x hF 2 ^[ 0 ,T] it follows that 

lim {J{v) - J(y^)) = 0. 


(3.86) 

(3.87) 
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Hence, we only need to prove that 


lim Xn{[v\n) = J{v) (3.88) 

n—>-oo 

The proof of (I3.88p is almost identical to the proof of Lemma 13.31 Lemma is proved. 

Having Lemmas 13.21 13.31 and 13.41 Theorem 11.21 follows from Lemma 12.21 
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